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Abstract: We investigate by means of a variational approach the effective potential of 
a I>D U{1) scalar model at finite temperature and compactified on and jZ 2 as well 
as the corresponding AD model obtained through a trivial dimensional reduction. We are 
particularly interested in the behavior of the thermal masses of the scalar field with respect 
to the Wilson line phase and the results obtained are compared with those coming from 
a one-loop effective potential calculation. We also explore the nature of the phase transition. 
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1 Introduction 


Studies on field theories at finite temperature have had a regain of interest particularly to 
investigate the electroweak transition and baryogenesis both in the Standard Model (SM) 
as effective low energy theory in 4 dimensions and in its extensions in five dimensions 
[1“5]. These new approaches have some of their theoretical roots in finite temperature 
held theory [6-9]. Considering topological aspects of a thermal formalism, it is realized 
that the prescription results in a scheme of compactihcation in time of the T = 0 theory. 
That is, the Matsubara results are equivalent to a path-integral calculated on x 

where is a circle of circumference (3 = 1/T. In the case of an extra dimension 
corresponding to 5'^/Z2 orbifold one should also take into account the parity projection 
[10]. As a consequence, the Matsubara prescription can be thought, in a generalized way, 
as a mechanism to deal simultaneously with spatial constraints and thermal effects in a 
held theory model. This concept has been developed for the Matsubara formalism by 
considering R^~'^ x x • • • x 5*"^, with li corresponding to inverse temperature and 
h corresponding to compactihcation of d — 1 spatial dimensions. Then the Feynman 

rules are modihed by introducing a generalized Matsubara prescription, performing the 
following multiple replacements (compactihcation of a d-dimensional subspace) [11], 
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i = 2, 3..., d — 1, 

( 1 . 1 ) 


where L,, i = 2, 3..., d — 1 are the sizes of the compactihed spatial dimensions. 

The ideas above have found new and interesting applications in the context of held 
theories with extra spatial dimensions, where the Higgs held is identihed as the zero mode 
of the hfth component of the gauge held, as for instance in [10, 12-27]. Models of this type 
are sometimes called models with gauge-Higgs unihcation. Perhaps these theories provide 
an interesting framework for physics beyond the SM. 

In particular, there has been a resurgence of interest in these ideas, as a new way to 
investigate the electroweak transition and baryogenesis. The electroweak phase transition, 
in 5-dimensional hnite-temperature held theory with a compactihed extra dimension, has 
been investigated in [2]. These authors compare theirs results with those obtained with the 
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SM at one-loop approximation. The conclusion is for a first-order transition with a strength 
inversely proportional to the Higgs mass, using the estimated value of the Higgs mass, of 
Mh < Mw/2. However, these authors also claim that when fermions in the uncompactified 
5-D space are introduced, more realistic values of the Higgs mass are obtained. In this case 
the first-order phase transition becomes weaker. Another interesting result of [2] is that 
up to temperatures of the order of T ~ 1/T, reliable (low order) perturbative calculations 
lead to reasonable results. It has been recently noticed that a 50(5) x 0(1) gauge Higgs 
unihcation in the Randall Sundrum metric can reproduce the Higgs mass at 125 GeV for 
three fermion spinorial representations and shown that the thermal phase transition at one 
loop is first order but very weak so baryogenesis could not occur [5]. 

From a physical and phenomenological point of view, an interest in theories with extra 
compactified dimensions at the inverse TeV scale arose in connection with the new LHC 
experiments. 

In order to go beyond one-loop approximation, self-consistent approaches have been 
also considered. One example is the Cornwall-Jackiw-Tomboulis effective action for com¬ 
posite operators that has been generalized to finite temperature [28-30]. Here we follow 
an alternative route, which consists in developing another, non perturbative, variational 
technique, the Gaussian Effective Potential (CEP) [31, 32], for both finite temperature and 
compactified spatial dimensions applying it as a tool to investigate effective field theories. 

The aim of this paper is to study, by using the CEP, the thermal masses of two models. 
One being the reduced model, at finite temperature, obtained from neglecting all mas¬ 
sive Kaluza Klein modes for a 5-dimensional U{1) model, thus trivially reducing it to a 
4-dimensional model. The other being the full 5-dimensional U{1) model, at finite temper¬ 
ature, compactified over 5^ and over the 5^/Z2 orbifold. Some of our results conhrm the 
perturbative results of [2]. 

In Section 2 we consider first the CEP for the reduced model where a truncation of the 
Kaluza Klein expansion has been performed, so that only the first KK mode is retained. 
In Section 3 we consider the full 5D scalar electrodynamics, deriving the CEP at finite 
temperature and calculating the Higgs thermal mass. We also discuss the structure of the 
phase transition by looking into the high temperature limit. 

2 The reduced model 

In this Section we review the reduced model which provides a simple exercise for studying 
the effective potential for the hve dimensional scalar electrodynamics by truncating the 
Kaluza Klein tower and considering only a first scalar lower mode. Let us consider the 
following lagrangian: 

£ = (2-1) 

which corresponds to the reduced 4Z? model. 

We introduce real fields ^i, 02, 


and perform a shift on the scalar field s, 


s(x) = s(x) + Sc 


(2.3) 


which allows to interpret Sc as a constant background field. 

In the canonically quantized version of the Gaussian effective potential (GEP) [31, 32] 
the quantum free fields cj) and s, respectively of masses and A, are expanded in the form, 


with 


and 
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a^{k)e-^^^ + h.c 


(2.4) 


(2.5) 


with ujf, = A:^ +A^. Standard commutation relations for creation and annihilation operators 
are assumed for os, 


[aj]{k),aj:{k')^ = 2ujk{27r)^5^{k - k') 


( 2 . 6 ) 


where S stands for either or A. In the Gaussian effective potential approach, the effective 
potential is evaluated as the expectation value of the Hamiltonian 


HG(a,f^,A) =< 0|^|0 > 


(2.7) 


where the vacuum |0) is annihilated by an{k) and a^{k) and % is the total Hamiltonian 
associated to the lagrangian density (2.1). Also we have redefined Sc = a/gR, 27rR being 
the length of the fifth compactified dimension . One can then minimize Vq with respect to 
the parameters S obtaining the associated gap equations which provide the values of the 
parameters that, when replaced in Vq, gives the Gaussian effective potential Vq 

An equivalent calculation may be performed by use of the first order (5-expansion [33, 
34], which amounts to split the Lagrangian C into a sum of a "free field" Cq term and a Cs 
term containing the interactions, such that when (5=1 the original Lagrangian is recovered. 
In the next section it will be convenient to cast our calculations into this approach. In both 
cases, the result for the reduced model is: 
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(2.9) 


with u!p = yp^ + 5 ^ Jq and R integrals are equivalent to the covariant form (in the 
Euclidean and for the R case up to an infinite constant) 
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The gap equations are: 


= ^ + XIom + g^Io{A), ( 2 . 11 ) 

= 2g^Ioin). 


where use has been done of the identity 

A/,(s) = S/o(S) (2.12) 

Replacing the values of fl and A from Eq.(2.11) into Eq.(2.8) gives the ’optimized’ 
Gaussian effective potential, Vq- The divergent integrals Iq,Ii are evaluated using dimen¬ 
sional regularization with the minimal subtraction scheme (MS) prescription. In particular 
one has (in the Euclidean and using their covariant form) 

= (^ ('“V -’) 


i = il? (‘°V -1) 

where /i is a regularization scale. Using this regularization prescription, Eq.(2.12) is still 
valid and therefore the form of the gap equation is the same as Eq.(2.11). 

The hnite temperature result can be obtained by replacing the I integrals by their finite 
temperature version as described in [35, 36]. For one has 

= + (2.15) 


where the parameter /3 indicates, as usual, the inverse of the temperature T and 
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(2.16) 

with c = 3/2 -|- 2 (log 47 r — 7 ) ~ 5.41. It is important to observe that the structure of the 

Gaussian effective potential is maintained as Eq.(2.12) is also valid for the hnite temperature 

integral With these substitutions one is then led to a hnite temperature Gaussian 
—FT 

effective potential Vq . 

Having obtained the optimized hnite temperature Gaussian effective potential we can 
calculate the Wilson line thermal mass of the held a, i.e.. 


M\T,a)=M-^, (2.17) 

where W is a suitable normalization factor. 

As we are dealing with an effective theory, a series of approximations are in order and 
so, following [2] we shall assume that g is negligible. In that case, we can consider that, 
from Eq.( 2 . 11 ), A = 0, and within the dimensional regularization scheme one can set scale- 
independent integrals as /o(0) or -^i(O) equal to zero. So we are led to a simpler situation 
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Figure 1. Squared thermal mass of the field a in the reduced model from the GEP (continuous 
line), one-loop effective potential (upper dashed line), improved one-loop (lower dashed line) for 
Lr= 1, A= 1, l/R. 


where only the equation for O needs to be satisfied. Assuming that LT = 1, A = 1 and 
H = l/R, the scale above which one considers that the effective theory breaks down, we 
plot M^{T, a), in Fig.l, for Vq (red line) and compare it with the results for the improved 
one-loop effective potential (red dashed line) and the "maive one-loop effective potential 
(upper dashed line) taken from [2], which, in our notation, is simply proportional to the 
integral /f'^(f2) for a value Al = a/R. The normalization factor Af is chosen such as to 
make equal to 1 the thermal mass calculated from the naive one-loop effective potential 
when a = 0. The main result is that the thermal mass calculated with GEP is substantially 
in agreement (slightly higher) with the prediction of the improved one-loop method. 


3 Compactified 5D scalar QED 


Let us now consider 5D U{1) scalar electrodynamics, with the 5th dimension compactified 
on a circle of length L = 2ttR, defined in a Euclidean space by the action 


p7ri? 


1 


Se = J dy j d'^x -EmatTmat + (T'm<p)^(T>m<p) 

+ + ^ 3 / 


with Cgf being the gauge hxing term and 


(3.1) 


Dm^ = + ig^AlMy^, 


(3.2) 


where tp satisfies the boundary condition (p{2ttR) = <^(0) and capital letters indicate 5D 
indices. Notice that the change from Minkowskian to Euclidean space corresponds to 
VmT^ _>■ —VmTm, therefore after a redefinition of sign in the action. We adopt the 
convention that capital Roman letter label indices running from 1 to 5, while Greek letters 
are used for indices which run from 1 to 4. We are interested in considering the effective 
potential for the vacuum solution for the 5th component of the field Am- Therefore we 
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will consider the Am held as a background field with only ^ 0 , where we have 

used the notation x = (xi, • • • 5 X 4 ) for the coordinates of the position vector in the non- 
compactified subspace. Here, as mentioned, it will be more convenient to introduce the 
Gaussian effective potential as the first order in the 5—expansion [33, 34], 

Introducing the real components of (p, 


if = ^(<^1 +f(^ 2 ), 


the part of the Lagrangian quadratic in can be rewritten as 

2 


with <I> = ( 01 , ^ 2 ) 


( ‘^gb^MdM \ 

-2g5AMdM -Om + gl^l j 


Then we shift the field A^ by a constant background field getting. 


(3.3) 


(3.4) 


(3.5) 


^5(x,X5) = H5(x,X5) - A 5 


(3.6) 


and we split the lagrangian as 

^ = (-^0 + ^int)s=l ( 3 - 7 ) 

where Cq is the sum of the quadratic lagrangian for (pip, the electromagnetic field terms 
and the terms containing the variational parameter 14: 


^0 — 




with 

p-i ^ ( ~^M + ^ + 2|(95 \ 

V “^ 1^5 -Qm + ^ +142 y ’ 

where we have used 55^5 = a/R as well as A 5 = LA. 

The interaction Lagrangian is then: 



(3.8) 


(3.9) 


(3.10) 


In Eq.(3.10), for simplicity, we kept only the quadratic terms that survive after the func¬ 
tional integrations. 

Introducing sources for the fields, the generating functional is written as: 


Z[J,j] = / D 01,02,-45 


exp 


/7 vR P 

<4X5 / 
-ttR J 


(Ax{Co+Cint} + (T 4 >) (j.H5(x,X5)) 


(3.11) 


6 








and the effective action is obtained by the Legendre transformation 


r[a] = In Z[j] - (j.A^), 


(3.12) 


where the brackets () are a shorthand notation for the integral f dx^ d^^x. 

One can expand the action, Eq.(3.12), to first order in 5 by the usual functional tech¬ 
niques getting, 


^GEP = ^Tr In {V ^)+^Tr In {G^ .^ ) > exp\^{< j'^.V.J > + < jG^j >)] 

(3.13) 

where, again, the brackets are a short notation for the integrations and the summations 
over the compactified variables, ) means that we have replaced the fields t^i, ^2 and 

A 5 by their functional derivatives in Ji, J 2 and j. The inverse of P in the momentum 
representation is given by 


Jd=o 


+ W + 2z> 5 y 

\ “2z^p5 Pm + / 


(3.14) 


The result, after the calculation of the traces and omitting the terms in g (as mentioned 
at the end of the previous section, g is considered negligible), reduces to 

Eg = hm - (3.15) 

where the integrals Iq and Ii, prior to compactification are given by: 
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The gap equation, obtained by minimizing the potential Vg with respect to the parameter 
P, in this case is: 

q,2 

= -^ + 2LXIo{n). (3.18) 

As before replacing the value of fl Eq.(3.18) into Eq.(3.15) gives the value of Gaussian 
effective potential, f^G) now for the compactified scalar QED. Formally, when the value Q 
gives a minimum, it is possible to replace Eq.(3.18) in Eq.(3.15), such that: 


VG = hm-—iiin). 


(3.19) 
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In order to introduce simultaneously the compactification of the 5th dimension over 
and a finite temperature, one has to perform the replacements, 
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Remembering that L = 27:R, we finally get. 
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In Eqs. (3.20) to (3.22) the sums over n and k refer respectively to the Matsubara and 
compactification modes. The motivation to go through this procedure, in order to obtain 
Eq.(3.22), is that this expression is much easier to manipulate by means of the Poisson 
summation formula. 
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together with the identity 
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Using Eqs. (3.23), (3.24) and the representation of Bessel functions of the third kind, 
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the integrals, and subtracting the divergent term corresponding to the zero mode, 
which however does not depend on a [2], can be written as: 
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In the calculations above we have performed the Poisson resummation over both the 
Matsubara [tj and compactification (k) modes. Equivalently, it is possible to use this 
procedure only over the Matsubara modes or only over the compactification ones. The 
results are: 
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(3.31) 


for the Matsubara modes. 

We have then laid the setup for the calculation of the thermal mass of the field a. In 
Fig. 2 we compare the results of the squared thermal mass obtained with the GEP (red 
continuous line), for LT = 1 and A = 1, with the improved one-loop effective potential 
[2]. The variational calculation result is larger than the improved one-loop for small values 
of a and smaller for larger value. In fact, the variational calculation interpolates the naif 
one-loop and improved results for the region of small values of a and is more sensitive than 
the perturbative approaches in the larger values of a region. 

It is also interesting to notice that, as in [2] we can obtain an indication of the order of 
the phase transition. Using the expressions of eqs.(3.28) and (3.29) for the end point 12 = 0 
and dropping all the modes, but the zero mode, in the Matsubara expansion we arrive at: 
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Figure 2. Squared thermal mass of the field a in the 5D model using the GEP (continuous line) 
and the improved one-loop (dashed line) for LT = 1 and A = 1. 


where we have also used 

The cosinus function in Eqs.(3.32) and (3.33) lead to ill-defined, rapidly oscillating 
series. Nevertheless, there is a solution in the framework of the zeta-function regulariza¬ 
tion. This procedure has already been used in [2] and it is a crucial step to these authors 
conclude for a first-order electroweak transition. As it is well-known, the zeta function 
can be analytically extended to the whole complex plane, having only one pole at z = 1. 
This analytical extension with a strictly negative even argument vanishes: (^{—2n) = 0 for 
integer n > 1. Accordingly, we expand the cosinus functions above in a power series of a 
and using = C(~2n) = 0 for all positive integers n > 1 it is possible to get, for 

definite values of q (we use ("(0) = —1/2)): 

-qa + {qaf) (3.34) 

and 

Now, the reasoning for doing these manipulations is that the = 0 endpoint corre¬ 
sponds to an infrared limit [2, 31], also one can see that in eqs.(3.28) and (3.29) when the 
temperature, T, is very large and for any non null i , the generalized Bessel function Bx 
approaches zero and the extremum values for the integrals occurs at the Q = 0 end point. 
Under these conditions the potential will be given by replacing eqs. (3.34) and (3.35) in 
Eq.(3.15) taking the endpoint = 0. 

In the compactification case we keep all degrees of freedom and the odd powers of 
a cancel out . One gets, apart from a term independent of ce, a polynomial of order 4 in a 
for the potential, which has the form, 

VG = ^^{Ac?^Ba^) (3.36) 

where the coefficients A = 33p -|- A/9^/(327r^) and B = 38p + 3A/9^/(327r^) are positive and 
p = LT. Since we are in the high temperature regime TL » 1, we expect the system to 
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be in the disordered phase. So, for A positive and when we keep all degrees of freedom, the 
form of Eq.(3.36) indicates the system has undergone a second-order phase transition. 

The IZ 2 compactification runs in a similar way, with the caveat that one has to 
project the states over a definite value for the charge q [10], we have to choose only one 
value of q for the expressions in Eqs. (3.34) and (3.35); in this case the odd powers of a 
do not cancel. That means that an term remains in the effective potential. Its form is 
proportional to 0^(35/48/9 -|- A/0^/(327r)^ indicating that for positive values of A the system 
has undergone a first order phase transition. 

Acknowledgments 

I.R. thanks the warm hospitality of the Dipartimento di Fisica Universita di Firenze where 
parts of this work were done. He also acknowledges fruitful discussions as well as cru¬ 
cial suggestions from Adolfo P.C. Malbouisson. The authors acknowledge partial financial 
support from INFN, Compaq and PRIN contract 2010YJ2NYW. 

References 

[1] C. Grojean, G. Servant and J. D. Wells, Phys. Rev. D71, 036001 (2005), [hep-ph/0407019]. 

[2] G. Panico and M. Serone, JHEP 05, 024 (2005), [hep-ph/0502255]. 

[3] N. Maru and K. Takenaga, Phys. Rev. D72, 046003 (2005), [hep-th/0505066]. 

[4] G. Delaunay, G. Grojean and J. D. Wells, JHEP 04, 029 (2008), [0711.2511]. 

[5] H. Hatanaka, 1304.5104. 

[6] T. Matsubara, Prog. Theor. Phys. 14, 351 (1955). 

[7] H. Ezawa, Y. Tomozawa and H. Umezawa, Nuovo Gim. 5, 810 (1957). 

[8] R. Kubo, J. Phys. Soc. Jap. 12, 570 (1957). 

[9] P. G. Martin and J. S. Schwinger, Phys. Rev. 115, 1342 (1959). 

[10] M. Kubo, G. Lim and H. Yamashita, Mod.Phys.Lett. A17, 2249 (2002), [hep-ph/0111327]. 

[11] F. G. Khanna, A. P. Malbouisson, J. M. Malbouisson and A. R. Santana, (2009), World 
Scientific, New Jersey, 2009 (ISBN-13: 978-981-281-887-4, ISBN-10: 981-281-887-1, ebook 
ISBN-13: 978-981-281-889-8, ebook ISBN-10: 981-281-889-8). 

[12] D. Fairlie, J.Phys.G G5, L55 (1979). 

[13] D. Fairlie, Phys.Lett. B82, 97 (1979). 

[14] N. Manton, Nucl.Phys. B158, 141 (1979). 

[15] Y. Hosotani, Phys.Lett. B126, 309 (1983). 

[16] Y. Hosotani, Phys.Lett. B129, 193 (1983). 

[17] 1. Antoniadis, Phys.Lett. B246, 377 (1990). 

[18] G.-L. Ho and Y. Hosotani, Nucl. Phys. B345, 445 (1990). 

[19] G. Dvali, S. Randjbar-Daemi and R. Tabbash, Phys.Rev. D65, 064021 (2002), 
[hep-ph/0102307]. 


11 


[20] L. J. Hall, Y. Nomura and D. Tucker-Smith, Nucl.Phys. B639, 307 (2002), [hep-ph/0107331]. 

]21] N. Arkani-Hamed, A. G. Cohen and H. Georgi, Phys.Lett. B516, 395 (2001), 
[hep-th/0103135]. 

]22] G. Burdman and Y. Nomura, Nucl.Phys. B656, 3 (2003), ]hep-ph/0210257]. 

]23] K. Agashe, R. Contino and A. Pomarol, Nucl.Phys. B719, 165 (2005), ]hep-ph/0412089]. 

]24] M. Serone, AIP Conf.Proc. 794, 139 (2005), [hep-ph/0508019]. 

[25] G. Panico, M. Serone and A. Wulzer, Nucl.Phys. B739, 186 (2006), ]hep-ph/0510373]. 

[26] Y. Hosotani, K. Oda, T. Ohnuma and Y. Sakamura, Phys.Rev. D78, 096002 (2008), 
[0806.0480]. 

[27] G. Panico, M. Safari and M. Serone, JHEP 1102, 103 (2011), [1012.2875]. 

[28] A. Barducci, R. Casalbuoni, D. Dominici, R. Gatto and G. Pettini, Phys. Lett. B179, 275 
(1986). 

[29] G. Amelino-Camelia and S.-Y. Pi, Phys. Rev. D47, 2356 (1993), [hep-ph/9211211[. 

[30] G. Smet, T. Vanzielighem, K. Van Acoleyen and H. Verschelde, Phys. Rev. D65, 045015 
(2002), [hep-th/0108163[. 

[31] P. M. Stevenson, Phys. Rev. D30, 1712 (1984). 

[32] P. M. Stevenson, Phys. Rev. D32, 1389 (1985). 

[33] A. Okopinska, Phys. Rev. D35, 1835 (1987). 

[34] I. Stancu and P. M. Stevenson, Phys. Rev. D42, 2710 (1990). 

[35] I. Roditi, Phys. Lett. B169, 264 (1986). 

[36] G. A. Hajj and P. M. Stevenson, Phys. Rev. D37, 413 (1988). 


12 - 


